arXiv:1501.01378v2 [hep-ex] 29 Apr 2015 


Measurement of yep in D (] — D [) oscillation using quantum correlations in 

e + e~ —>■ D°D° at yfs = 3.773 GeV 


M. Ablikim 1 , M. N. Achasov 8 ’ 0 , X. C. Ai 1 , O. Albayrak 4 , M. Albrecht 3 , D. J. Ambrose 43 , A. Amoroso 47A ’ 47 °, 

F. F. An 1 , Q. An 44 , J. Z. Bai 1 , R. Baldini Ferroli 19A , Y. Ban 30 , D. W. Bennett 18 , J. V. Bennett 4 , M. Bertani 19A , 
D. Bettoni 20A , J. M. Bian 42 , F. Bianchi 47A ’ 47 °, E. Boger 22 ’\ O. Bondarenko 24 , I. Boyko 22 , R. A. Briere 4 , H. Cai 49 , 

X. Cai 1 , O. Cakir 39A ’ & , A. Calcaterra 19A , G. F. Cao 1 , S. A. Cetin 39S , J. F. Chang 1 , G. Chelkov 22 ’ c , G. Chen 1 , 
H. S. Chen 1 , H. Y. Chen 2 , J. C. Chen 1 , M. L. Chen 1 , S. J. Chen 28 , X. Chen 1 , X. R. Chen 25 , Y. B. Chen 1 , 

H. P. Cheng 16 , X. K. Chu 30 , G. Cibinetto 20A , D. Cronin-Hennessy 42 , H. L. Dai 1 , J. P. Dai 33 , A. Dbeyssi 13 , 

D. Dedovich 22 , Z. Y. Deng 1 , A. Denig 21 , I. Denysenko 22 , M. Destefanis 47A,47 °, F. De Mori 47A ’ 47 °, Y. Ding 26 , 

C. Dong 29 , J. Dong 1 , L. Y. Dong 1 , M. Y. Dong 1 , S. X. Du 51 , P. F. Duan 1 , J. Z. Fan 38 , J. Fan g 1 , S. S. Fang 1 , 

X. Fang 44 , Y. Fang 1 , L. Fava 47S ’ 47 °, F. Feldbauer 21 , G. Felici 19A , C. Q. Feng 44 , E. Fioravanti 20A , M. Fritsch 13 ’ 21 , 

C. D. Fu 1 , Q. Gao 1 , Y. Gao 38 , Z. Gao 44 , I. Garzia 20A , K. Goetzen 9 , W. X. Gong 1 , W. Gradl 21 , M. Greco 47A ’ 47 °, 
M. H. Gu 1 , Y. T. Gu 11 , Y. H. Guan 1 , A. Q. Guo 1 , L. B. Guo 27 , T. Guo 27 , Y. Guo 1 , Y. P. Guo 21 , Z. Haddadi 24 , 
A. Hafner 21 , S. Han 49 , Y. L. Han 1 , F. A. Harris 41 , K. L. He 1 , Z. Y. He 29 , T. Held 3 , Y. K. Heng 1 , Z. L. Hou 1 , 

C. Hu 27 , H. M. Hu 1 , J. F. Hu 47A , T. Hu 1 , Y. Hu 1 , G. M. Huang 5 , G. S. Huang 44 , H. P. Huang 49 , J. S. Huang 14 , 
X. T. Huang 32 , Y. Huang 28 , T. Hussain 46 , Q. Ji 1 , Q. P. Ji 29 , X. B. Ji 1 , X. L. Ji 1 , L. L. Jiang 1 , L. W. Jiang 49 , 

X. S. Jiang 1 , J. B. Jiao 32 , Z. Jiao 16 , D. P. Jin 1 , S. Jin 1 , T. Johansson 48 , A. Julin 42 , N. Kalantar-Nayestanaki 24 , 
X. L. Kang 1 , X. S. Kang 29 , M. Kavatsyuk 24 , B. C. Ke 4 , R. Kliemt 13 , B. Kloss 21 , O. B. Kolcu 39S ’ d , B. Kopf 3 , 

M. Kornicer 41 , W. Kuehn 23 , A. Kupsc 48 , W. Lai 1 , J. S. Lange 23 , M. Lara 18 , P. Larin 13 , C. H. Li 1 , Cheng Li 44 , 
D. M. Li 51 , F. Li 1 , G. Li 1 , H. B. Li 1 , J. C. Li 1 , Jin Li 31 , K. Li 32 , K. Li 12 , P. R. Li 40 , T. Li 32 , W. D. Li 1 , W. G. Li 1 , 

X. L. Li 32 , X. M. Li 11 , X. N. Li 1 , X. Q. Li 29 , Z. B. Li 37 , H. Liang 44 , Y. F. Liang 35 , Y. T. Liang 23 , G. R. Liao 10 , 

D. X. Lin 13 , B. J. Liu 1 , C. X. Liu 1 , F. H. Liu 34 , Fang Liu 1 , Feng Liu 5 , H. B. Liu 11 , H. H. Liu 1 , H. H. Liu 15 , 

H. M. Liu 1 , J. Liu 1 , J. P. Liu 49 , J. Y. Liu 1 , K. Liu 38 , K. Y. Liu 26 , L. D. Liu 30 , P. L. Liu 1 , Q. Liu 40 , S. B. Liu 44 , 

X. Liu 25 , X. X. Liu 40 , Y. B. Liu 29 , Z. A. Liu 1 , Zhiqiang Liu 1 , Zhiqing Liu 21 , H. Loehner 24 , X. C. Lou 1 ’ 0 , 

H. J. Lu 16 , J. G. Lu 1 , R. Q. Lu 17 , Y. Lu 1 , Y. P. Lu\ C. L. Luo 27 , M. X. Luo 50 , T. Luo 41 , X. L. Luo 1 , M. Lv\ 

X. R. Lyu 40 , F. C. Ma 26 , H. L. Ma 1 , L. L. Ma 32 , Q. M. Ma 1 , S. Ma 1 , T. Ma 1 , X. N. Ma 29 , X. Y. Ma 1 , 

F. E. Maas 13 , M. Maggiora 47A ’ 47 °, Q. A. Malik 46 , Y. J. Mao 30 , Z. P. Mao 1 , S. Marcello 4721 ’ 470 , 

J. G. Messchendorp 24 , J. Min 1 , T. J. Min 1 , R. E. Mitchell 18 , X. H. Mo 1 , Y. J. Mo 5 , C. Morales Morales 13 , 

K. Moriya 18 , N. Yu. Muchnoi 8 ’ 0 , H. Muramatsu 42 , Y. Nefedov 22 , F. Nerling 13 , I. B. Nikolaev 8 ’ 0 , Z. Ning 1 , 

S. Nisar 7 , S. L. Niu 1 , X. Y. Niu 1 , S. L. Olsen 31 , Q. Ouyang 1 , S. Pacetti 19S , P. Patteri 19A , M. Pelizaeus 3 , 

H. P. Peng 44 , K. Peters 9 , J. L. Ping 27 , R. G. Ping 1 , R. Poling 42 , Y. N. Pu 17 , M. Qi 28 , S. Qian 1 , C. F. Qiao 40 , 

L. Q. Qin 32 , N. Qin 49 , X. S. Qin 1 , Y. Qin 30 , Z. H. Qin 1 , J. F. Qiu 1 , K. H. Rashid 46 , C. F. Redmer 21 , H. L. Ren 17 , 

M. Ripka 21 , G. Rong 1 , X. D. Ruan 11 , V. Santoro 20A , A. Sarantsev 22 ’-^, M. Savrie 20B , K. Schoenning 48 , 

S. Schumann 21 , W. Shan 30 , M. Shao 44 , C. P. Shen 2 , P. X. Shen 29 , X. Y. Shen 1 , H. Y. Sheng 1 , M. R. Shepherd 18 , 
W. M. Song 1 , X. Y. Song 1 , S. Sosio 47A ’ 47 °, S. Spataro 47A ’ 47 °, B. Spruck 23 , G. X. Sun 1 , J. F. Sun 14 , S. S. Sun 1 , 

Y. .J. Sun 44 , Y. Z. Sun 1 , Z. .J. Sun 1 , Z. T. Sun 18 , C. J. Tang 35 , X. Tang 1 , I. Tapan 390 , E. H. Thorndike 43 , 

M. Tiemens 24 , D. Toth 42 , M. Ullrich 23 , I. Uman 39B , G. S. Varner 41 , B. Wang 29 , B. L. Wang 40 , D. Wang 30 , 

D. Y. Wang 30 , K. Wang 1 , L. L. Wang 1 , L. S. Wang 1 , M. Wang 32 , P. Wang 1 , P. L. Wang 1 , Q. J. Wang 1 , 

S. G. Wang 30 , W. Wang 1 , X. F. Wang 38 , Y. D. Wang 19A , Y. F. Wang 1 , Y. Q. Wang 21 , Z. Wang 1 , Z. G. Wang 1 , 

Z. H. Wang 44 , Z. Y. Wang 1 , T. Weber 21 , D. H. Wei 10 , J. B. Wei 30 , P. Weidenkaff 21 , S. P. Wen 1 , U. Wiedner 3 , 

M. Wolke 48 , L. H. Wu\ Z. Wu 1 , L. G. Xia 38 , Y. Xia 17 , D. Xiao 1 , Z. J. Xiao 27 , Y. G. Xie 1 , G. F. Xu 1 , L. Xu 1 , 

Q. J. Xu 12 , Q. N. Xu 40 , X. P. Xu 36 , L. Yan 44 , W. B. Yan 44 , W. C. Yan 44 , Y. H. Yan 17 , H. X. Yang 1 , L. Yang 49 , 

Y. Yang 5 , Y. X. Yang 10 , H. Ye 1 , M. Ye 1 , M. H. Ye 6 , J. H. Yin 1 , B. X. Yu 1 , C. X. Yu 29 , H. W. Yu 30 , J. S. Yu 25 , 
C. Z. Yuan 1 , W. L. Yuan 28 , Y. Yuan 1 , A. Yuncu 395 ’ 9 , A. A. Zafar 46 , A. Zallo 19A , Y. Zeng 17 , B. X. Zhang 1 , 

B. Y. Zhang 1 , C. Zhang 28 , C. C. Zhang 1 , D. H. Zhang 1 , H. H. Zhang 37 , H. Y. Zhang 1 , J. J. Zhang 1 , J. L. Zhang 1 , 
J. Q. Zhang 1 , J. W. Zhang 1 , J. Y. Zhang 1 , J. Z. Zhang 1 , K. Zhang 1 , L. Zhang 1 , S. H. Zhang 1 , X. Y. Zhang 32 , 

Y. Zhang 1 , Y. H. Zhang 1 , Y. T. Zhang 44 , Z. H. Zhang 5 , Z. P. Zhang 44 , Z. Y. Zhang 49 , G. Zhao 1 , J. W. Zhao 1 , 

J. Y. Zhao 1 , J. Z. Zhao 1 , Lei Zhao 44 , Ling Zhao 1 , M. G. Zhao 29 , Q. Zhao 1 , Q. W. Zhao 1 , S. J. Zhao 51 , 

T. C. Zhao 1 , Y. B. Zhao 1 , Z. G. Zhao 44 , A. Zhemchugov 22 ’' 1 , B. Zheng 45 , J. P. Zheng 1 , W. J. Zheng 32 , 

Y. H. Zheng 40 , B. Zhong 27 , L. Zhou 1 , Li Zhou 29 , X. Zhou 49 , X. K. Zhou 44 , X. R. Zhou 44 , X. Y. Zhou 1 , K. Zhu 1 , 
K. J. Zhu 1 , S. Zhu 1 , X. L. Zhu 38 , Y. C. Zhu 44 , Y. S. Zhu 1 , Z. A. Zhu 1 , J. Zhuang 1 , B. S. Zou 1 , J. H. Zou 1 

(BESIII Collaboration) 


1 


1 Institute of High Energy Physics, Beijing 100049, People’s Republic of China 
2 Beihang University, Beijing 100191, People’s Republic of China 
3 Bochum Ruhr-University, 0-44180 Bochum, Germany 
4 Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA 
5 Central China Normal University, Wuhan 430079, People’s Republic of China 
6 China Center of Advanced Science and Technology, Beijing 100190, People’s Republic of China 
7 COMSATS Institute of Information Technology, Lahore, Defence Road, Off Raiwind Road, 54000 Lahore, 

Pakistan 

8 G.I. Budker Institute of Nuclear Physics SB RAS (BINP), Novosibirsk 630090, Russia 
9 GSI Helmholtzcentre for Heavy Ion Research GmbH, D-64291 Darmstadt, Germany 
10 Guangxi Normal University, Guilin 541004, People’s Republic of China 
11 GuangXi University, Nanning 530004, People’s Republic of China 
12 Hangzhou Normal University, Hangzhou 310036, People’s Republic of China 
13 Helmholtz Institute Mainz, Johann-Joachim-Becher-Weg 45, D-55099 Mainz, Germany 
14 Henan Normal University, Xinxiang 453007, People’s Republic of China 
15 Henan University of Science and Technology, Luoyang 471003, People’s Republic of China 
16 Huangshan College, Huangshan 245000, People’s Republic of China 
17 Hunan University, Changsha 410082, People’s Republic of China 
18 Indiana University, Bloomington, Indiana 47405, USA 
19 (A)INFN Laboratori Nazionali di Frascati, 1-00044> Frascati, Italy; (B)INFN and University of Perugia, 

1-06100, Perugia, Italy 

20 (A)INFN Sezione di Ferrara, 1-441%%, Ferrara, Italy; (B)University of Ferrara, I~441%%, Ferrara, Italy 

21 Johannes Gutenberg University of Mainz, Johann-Joachim-Becher-Weg 45, D-55099 Mainz, Germany 

22 Joint Institute for Nuclear Research, 141980 Dubna, Moscow region, Russia 
23 Justus Liebig University Giessen, II. Physikalisches Institut, Heinrich-Buff-Ring 16, D-35392 Giessen, Germany 
24 KVI-CART, University of Groningen, NL-9747 AA Groningen, The Netherlands 
25 Lanzhou University, Lanzhou 730000, People’s Republic of China 
26 Liaoning University, Shenyang 110036, People’s Republic of China 
27 Nanjing Normal University, Nanjing 210023, People’s Republic of China 

28 Nanjing University, Nanjing 210093, People’s Republic of China 

29 Nankai University, Tianjin 300071, People’s Republic of China 

30 Peking University, Beijing 100871, People’s Republic of China 

31 Seoul National University, Seoul, 151-747 Korea 
32 Shandong University, Jinan 250100, People’s Republic of China 
33 Shanghai Jiao Tong University, Shanghai 200240, People’s Republic of China 
34 Shanxi University, Taiyuan 030006, People’s Republic of China 

35 Sichuan University, Chengdu 610064, People’s Republic of China 

36 Soochow University, Suzhou 215006, People’s Republic of China 

37 Sun Yat-Sen University, Guangzhou 510275, People’s Republic of China 
38 Tsinghua University, Beijing 100084, People’s Republic of China 
39 (A)Istanbul Aydin University, 34%95 Sefakoy, Istanbul, Turkey; (B)Dogus University, 34722 Istanbul, Turkey; 

(C)Uludag University, 16059 Bursa, Turkey 

40 University of Chinese Academy of Sciences, Beijing 100049, People’s Republic of China 
41 University of Hawaii, Honolulu, Hawaii 96822, USA 
42 University of Minnesota, Minneapolis, Minnesota 55455, USA 
43 University of Rochester, Rochester, New York 14627, USA 
44 University of Science and Technology of China, Hefei 230026, People’s Republic of China 
45 University of South China, Hengyang 421001, People’s Republic of China 
46 University of the Punjab, Lahore-54590, Pakistan 

47 (A)University of Turin, 1-10125, Turin, Italy; (B)University of Eastern Piedmont, 1-15121, Alessandria, Italy; 

(C)INFN, 1-10125, Turin, Italy 
48 Uppsala University, Box 516, SE-75120 Uppsala, Sweden 
49 Wuhan University, Wuhan 430072, People’s Republic of China 
50 Zhejiang University, Hangzhou 310027, People’s Republic of China 
51 Zhengzhou University, Zhengzhou 450001, People’s Republic of China 

a Also at the Novosibirsk State University, Novosibirsk, 630090, Russia 


2 


b Also at Ankara University, 06100 Tandogan, Ankara, Turkey 
c Also at the Moscow Institute of Physics and Technology, Moscow 141700, Russia and at the Functional 
Electronics Laboratory, Tomsk State University, Tomsk, 634050, Russia 
d Currently at Istanbul Arel University, 34295 Istanbul, Turkey 
e Also at University of Texas at Dallas, Richardson, Texas 75083, USA 
f Also at the PNPI, Gatchina 188300, Russia 
9 Also at Bogazici University, 34342 Istanbul, Turkey 
h Also at the Moscow Institute of Physics and Technology, Moscow 141700, Russia 


Abstract 

We report a measurement of the parameter yep in D° — D° oscillations performed by taking advantage of 
quantum coherence between pairs of D°D° mesons produced in e + e _ annihilations near threshold. In this 
work, doubly-tagged D°D a events, where one D decays to a CP eigenstate and the other D decays in a 
semileptonic mode, are reconstructed using a data sample of 2.92 fb -1 collected with the BESIII detector 
at the center-of-mass energy of y/s = 3.773 GeV. We obtain yep = (—2.0 ± 1.3 ± 0.7)%, where the first 
uncertainty is statistical and the second is systematic. This result is compatible with the current world 
average. 

Keywords: BESIII, D° — D° oscillation, yep , quantum correlation 


1. Introduction 

1.1. Charm oscillation 

It is well known that oscillations between meson and antimeson, also called mixing, can occur when the 
flavor eigenstates differ from the physical mass eigenstates. These effects provide a mechanism whereby 
interference in the transition amplitudes of mesons and antimesons may occur. They may also allow for 
manifestation of CP violation ( CPV ) in the underlying dynamics [1, 2]. Oscillations in the K° — K° [3], 
B° — B° [4] and B® — B ° [5] systems are established; their oscillation rates are well-measured and consistent 
with predictions from the standard model (SM) [6]. After an accumulation of strong evidence from a variety 
of experiments [7, 8, 9], D° — D° oscillations were recently firmly established by LHCb [10]. The results 
were soon confirmed by CDF [11] and Belle [12]. 

The oscillations are conventionally characterized by two dimensionless parameters x = A m/T and y = 
AT/2T, where Am and AT are the mass and width differences between the two mass eigenstates and T is 
the average decay width of those eigenstates. The mass eigenstates can be written as |Di j2 ) = p\D°)±q\D Q ), 
where p and q are complex parameters and <f> = arg (q/p) is a CP-violating phase. Using the phase convention 
CP|.D 0 ) = +|D°), the CP eigenstates of the D meson can be written as 

lDcpp) JKl±fl, (1 , 

The difference in the effective lifetime between D decays to CP eigenstates and flavor eigenstates can be 
parameterized by yep- In the absence of direct CPV, but allowing for small indirect CPV, we have [13] 
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Table 1: D final states reconstructed in this analysis. 


Type 

Mode 

CP+ 

K + Ii~, 7r + 7r - , Ar§7r°7r° 

CP- 

KgTT 0 , AigW, KgT] 

Semileptonic 

K^e ± v 1 


In the absence of CPV 7 one has \p/q\ = 1 and <j> = 0, leading to yep = y- 

Although D° — D° mixing from short-distance physics is suppressed by the CKM matrix [14, 15] and the 
GIM mechanism [16], sizeable charm mixing can arise from long-distance processes and new physics [1, 17]. 
Current experimental precision [18] is not sufficient to conclude whether physics beyond the SM is involved, 
and further constraints are needed. So far, the most precise determination of the size of the mixing has been 
obtained by measuring the time-dependent decay rate in the D —> K ± tt t channel [10, 11, 12]. However, the 
resulting information on the mixing parameters x and y is highly correlated. It is important to access the 
mixing parameters x and y directly to provide complementary constraints. 

In this analysis, we use a time-integrated method to extract yep , as proposed in the references [19, 20, 
21, 22], which uses threshold D°D° pair production in e + e _ — > 7 * — > D°D°. In this process, the D°D° 
pair is in a state of definite C = —1, such that the two D mesons necessarily have opposite CP eigenvalues. 
The method utilizes the semileptonic decays of D meson and hence, avoids the complications from hadronic 
effects in D decays, thus provides a clean and unique way to probe the D a — D° oscillation. 


1.2. Formalism 


In the semileptonic decays of neutral D mesons (denoted as D —> l) 1 , the partial decay width is only 
sensitive to flavor content and does not depend on the CP eigenvalue of the parent D meson. However, 
the total decay width of the Dcp± does depend on its CP eigenvalue: Tqp± = T(1 ± yep )■ Thus, the 
semileptonic branching fraction of the CP eigenstates Dqp± is Bd cp± ^i ~ £(d->z(1 ~F Ugp ), and yep can 
be obtained as 


VCP 


1 /B DcP-^rl 
4 \Bp) cp+ ^l 


^Dcpp^l 

Bd cp _^i 


( 3 ) 


At BESIII, quantum-correlated D°D° pairs produced at threshold allow us to measure Bd cp± ^i- Specif¬ 
ically, we begin with a fully reconstructed D candidate decaying into a CP eigenstate, the so-called Single 
Tag (ST). We have thus tagged the CP eigenvalue of the partner D meson. For a subset of the ST events, 
the so-called Double Tag (DT) events, this tagged partner D meson is also observed via one of the semilep¬ 
tonic decay channels. CP violation in D decays is known to be very small [18], and can be safely neglected. 
Therefore, Bd cp ^.->i can be obtained as 


^Dcpzf —>Z 


Ncp±d £CP± 
Ncp± £cp±-.i ’ 


( 4 ) 


where Ncp± ( Ncp±-i ) and ecp± (ecp±-j) denote the signal yields and detection efficiencies of ST decays 
D — > CPF (DT decays DD —> CP±: l ), respectively. For CP eigenstates, as listed in Table 1, we choose 
modes with unambiguous CP content and copious yields. The CP violation in Kg decays is known to be 
very small, it is therefore neglected. The semileptonic modes used for the DT in this analysis are 
and RTyFv. 


1.3. The BESIII detector and data sample 

The analysis presented in this paper is based on a data sample with an integrated luminosity of 
2.92 fb -1 [23] collected with the BESIII detector [24] at the center-of-mass energy of y/s = 3.773 GeV. 


1 Charge-conjugate modes are implied. 
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The BESIII detector is a general-purpose solenoidal detector at the BEPCII [25] double storage rings. The 
detector has a geometrical acceptance of 93% of the full solid angle. We briefly describe the components of 
BESIII from the interaction point (IP) outwards. A small-cell main drift chamber (MDC), using a helium- 
based gas to measure momenta and specific ionizations of charged particles, is surrounded by a time-of-flight 
(TOF) system based on plastic scintillators that determines the flight times of charged particles. A CsI(Tl) 
electromagnetic calorimeter (EMC) detects electromagnetic showers. These components are all situated in¬ 
side a superconducting solenoid magnet, that provides a 1.0 T magnetic field parallel to the beam direction. 
Finally, a multi-layer resistive plate counter system installed in the iron flux return yoke of the magnet is 
used to track muons. The momentum resolution for charged tracks in the MDC is 0.5% for a transverse mo¬ 
mentum of 1 GeV/c. The energy resolution for showers in the EMC is 2.5% (5.0%) for 1 GeV photons in the 
barrel (end cap) region. More details on the features and capabilities of BESIII can be found elsewhere [24] . 

High-statistics Monte Carlo (MC) simulations are used to evaluate the detection efficiency and to un¬ 
derstand backgrounds. The GEANT4-based [26] MC simulation program is designed to simulate interactions 
of particles in the spectrometer and the detector response. For the production of *0(3770), the KKMC [27] 
package is used; the beam energy spread and the effects of initial-state radiation (ISR) are included. The 
MC samples consist of the DD pairs with consideration of quantum coherence for all modes relevant to this 
analysis, non -DD decays of 0(3770), ISR production of low-mass ip states, and QED and qq continuum 
processes. The effective luminosity of the MC samples is about 10 times that of the analyzed data. Known 
decays recorded by the Particle Data Group (PDG) [ 6 ] are generated with EVTGEN [28, 29] using PDG 
branching fractions, and the remaining unknown decays are generated with LUNDCHARM [30]. Final-state 
radiation (FSR) of charged tracks is taken into account with the PHOTOS package [31]. 

2. Event selection and data analysis 

Each charged track is required to satisfy | cos0| < 0.93, where 9 is the polar angle with respect to the 
beam axis. Charged tracks other than Kg daughters are required to be within 1 cm of the IP transverse to 
the beam line and within 10 cm of the IP along the beam axis. Particle identification for charged hadrons h 
( h = 7 r, K) is accomplished by combining the measured energy loss {dE/dx) in the MDC and the flight time 
obtained from the TOF to form a likelihood C(h) for each hadron hypothesis. The K± (tt^) candidates are 
required to satisfy C(K) > C(n) (£( tt) > C(K)). 

The Kg candidates are selected with a vertex-constrained fit from pairs of oppositely charged tracks, 
which are required to be within 20 cm of the IP along the beam direction; no constraint in the transverse 
plane is required. The two charged tracks are not subjected to the particle identification discussed above, 
and are assumed to be pions. We impose 0.487 GeV/c 2 < M n + n - < 0.511 GeV/c 2 , that is within about 
3 standard deviations of the observed Kg mass, and the two tracks are constrained to originate from a 
common decay vertex by requiring the y 2 of the vertex fit to be less than 100. The decay vertex is required 
to be separated from the IP with a significance greater than two standard deviations. 

Reconstructed EMC showers that are separated from the extrapolated positions of any charged tracks 
by more than 10 standard deviations are taken as photon candidates. The energy deposited in nearby TOF 
counters is included to improve the reconstruction efficiency and energy resolution. Photon candidates must 
have a minimum energy of 25 MeV for barrel showers (| cos 9\ < 0.80) and 50 MeV for end cap showers (0.84 < 

| cos0| < 0.92). The showers in the gap between the barrel and the end cap regions are poorly reconstructed 
and thus excluded. The shower timing is required to be no later than 700 ns after the reconstructed event 
start time to suppress electronic noise and energy deposits unrelated to the event. The r] and 7 r° candidates 
are reconstructed from pairs of photons. Due to the poorer resolution in the EMC end cap regions, those 
candidates with both photons coming from EMC end caps are rejected. The invariant mass M 77 is required 
to be 0.115 GeV/c 2 < M 77 < 0.150 GeV/c 2 for 7r° and 0.505 GeV/c 2 < M 77 < 0.570 GeV/c 2 for r/ 
candidates. The photon pair is kinematically constrained to the nominal mass of the 7 r° or rj [ 6 ] to improve 
the meson four-vector calculation. 

The ui candidates are reconstructed through the decay ui —> 7 r + 7 r _ 7 r°. For all modes with u> candidates, 
sideband events in the M w + n - n o spectrum are used to estimate peaking backgrounds from non-w D —> 
KgTT + TT~TT° decays. We take the signal region as (0.7600, 0.8050) GeV/c 2 and the sideband regions as 


5 



Figure 1: Fit to the invariant mass for events reconstructed from data. The solid line is the total fit and the dashed 

line shows the polynomial background. The shaded area shows the signal region and cross-hatched areas show the sidebands. 


(0.6000, 0.7300) GeV/c 2 or (0.8300, 0.8525) GeV/c 2 . The upper edge of the right sideband is restricted 
because of the K*p background that alters the shape of The sidebands are scaled to the estimated 

peaking backgrounds in the signal region. The scaling factor is determined from a fit to the 
distribution in data, as shown in Fig. 1, where the u> signal is determined with the MC shape convoluted 
with a Gaussian whose parameters are left free in the fit to better match data resolution, and the background 
is modeled by a polynomial function. 

2.1. Single tags using CP modes 

To identify the reconstructed D candidates, we use two variables, the beam-constrained mass Mbc and 
the energy difference AE, which are defined as 

M BC = ^beam/c 4 - \Pd\ 2 /(?, (5) 

AE = Ed — -Ebeam, (6) 

where pd and Ed are the momentum and energy of the D candidate in the e + e _ center-of-mass system, 
and F^beam is the beam energy. The D signal peaks at the nominal D mass in Mbc and at zero in AE. 
We accept only one candidate per mode per event; when multiple candidates are present, the one with the 
smallest |Ai?| is chosen. Since the correlation between AE and Mbc is found to be small, this will not bias 
the background distribution in Mbc- We apply the mode-dependent AE requirements listed in Table 2. 

For K + K~ and 7r + 7r“ ST modes, if candidate events contain only two charged tracks, the following 
requirements are applied to suppress backgrounds from cosmic rays and Bhabha events. First, we require 
at least one EMC shower separated from the tracks of the ST with energy larger than 50 MeV. Second, the 
two ST tracks must not be both identified as muons or electrons, and, if they have valid TOF times, the 
time difference must be less than 5 ns. Based on MC studies, no peaking background is present in Mbc 
in our ST modes except for the Kgir 0 mode. In the Kg 7r° ST mode, there are few background events 
from D° —> pn. From MC studies, the estimated fraction is less than 0.5%; this will be considered in the 
systematic uncertainties. 

The Mbc distributions for the six ST modes are shown in Fig. 2. Unbinned maximum likelihood fits are 
performed to obtain the numbers of ST yields except in the Kgio mode, for which a binned least-square fit is 
applied to the Mbc distribution after subtraction of the w sidebands. In each fit, the signal shape is derived 
from simulated signal events convoluted with a bifurcated Gaussian with free parameters to account for 
imperfect modeling of the detector resolution and beam energy calibration. Backgrounds are described by the 
ARGUS [32] function. The measured ST yields in the signal region of 1.855 GeV/c 2 < Mbc <1-875 GeV/c 2 
and the corresponding efficiencies are given in Table 3. 
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Table 2: Requirements on AS for ST D candidates. 


Mode Requirement (GeV) 


K+K~ 

-0.020 

< 

A E 

< 

0.020 

7r + 7T _ 

-0.030 

< 

A E 

< 

0.030 

K%TT°n° 

-0.080 

< 

A E 

< 

0.045 

K° s n° 

-0.070 

< 

A E 

< 

0.040 


-0.050 

< 

A E 

< 

0.030 

K sV 

-0.040 

< 

A E 

< 

0.040 



Figure 2: The Mbc distributions for ST D candidates from data. The solid line is the total fit and the dashed line shows the 
background contribution described by an ARGUS function. 
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Table 3: Yields and efficiencies of all ST and DT modes, where Nqp± ( Ncp±i ) and £cp± ( £ CP±l ) denote signal yields and 
detection efficiencies of D —> CP zb (DD —> CP=t; /), respectively. The uncertainties are statistical only. 


ST Mode 

Nqp± 

£cp± (%) 

K+K~ 

54494 ± 251 

61.32 ± 0.18 

=1 

+ 

si 

i 

19921 ± 174 

64.09 ± 0.18 

Kgir°Tr° 

24015 ± 236 

16.13 ± 0.08 

K° s n° 

71421 ± 285 

40.67 ± 0.14 


20989 ± 243 

13.44 ± 0.07 


9878 ± 117 

34.39 ± 0.13 

DT Mode 

NcP±;l 

£CP±-,i (%) 

K+K-, Kev 

1216 ± 40 

39.80 ± 0.14 

7r + 7r _ , Kev 

427 ± 23 

41.75 ± 0.14 

Kgir°ir 0 , Kev 

560 ± 28 

11.05 ± 0.07 

Kgir 0 , Kev 

1699 ± 47 

26.70 ± 0.12 

KgUj, Kev 

481 ± 30 

9.27 ± 0.07 

Kgij, Kev 

243 ± 17 

22.96 ± 0.11 

K+K-, Kfjv 

1093 ± 37 

36.89 ± 0.14 

7r + 7r _ , K^lv 

400 ± 23 

38.43 ± 0.15 

K%TT°n°, Kfiv 

558 ± 28 

10.76 ± 0.08 

Kgir 0 , K[iv 

1475 ± 43 

25.21 ± 0.12 

KgUJ, Kfj,v 

521 ± 27 

8.75 ± 0.07 


241 ± 18 

21.85 ± 0.11 


2.2. Double tags of semileptonic modes 

In each ST event, we search among the unused tracks and showers for semileptonic D —> Ke(f£)v candi¬ 
dates. We require that there be exactly two oppositely-charged tracks that satisfy the fiducial requirements 
described above. 

In searching for K^jlv decays, kaon candidates are required to satisfy C(K) > £( 71 -). If the two tracks can 
pass the criterion, the track with larger £(/\) is taken as the K ± candidate, and the other track is assumed 
to be the /x candidate. The energy deposit in the EMC of the fi candidate is required to be less than 
0.3 GeV. We further require the K^i invariant mass to be less than 1.65 GeV/c 2 to reject D —> Kir 
backgrounds. The total energy of remaining unmatched EMC showers, denoted as £ ex t ra , is required to be 
less than 0.2 GeV to suppress D Kirn 0 backgrounds. To reduce backgrounds from the D Kev process, 
the ratio TZ C > (e) = £ (e)/[£ (e) + £ (/x) +£ (it) + £ (K)} is required to be less than 0.8, where the likelihood 
£ ( i ) for the hypothesis i = e, /x, it or K is formed by combining EMC information with the dE/dx and 
TOF information. 

To select Kev events, electron candidates are required to satisfy £ (e) > 0.001 and 7 Z c , (e) >0.8, where 

lZ c ,(e) = £ (e)/[£ (e)+£ (ir)+£ (K)]. If both tracks satisfy these requirements, the one with larger 1Z C , (e) 
is taken as the electron. The remaining track is required to satisfy £(K) > £(ir). 

The variable U m i ss is used to distinguish semileptonic signal events from background: 

bifiiss = E m ‘ lss c|prniss|, (7) 


where, 


Emiss — -Gbeam K E[, 


P miss — 



( 8 ) 

(9) 
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E K (i) ( Pk(i )) is the energy (three-momentum) of K T (l±), Pst is the unit vector in the direction of the 
reconstructed CP-tagged D and mu is the nominal D mass. The use of the beam energy and the nominal 
D mass for the magnitude of the CP-tagged D improves the P m i ss resolution. Since E equals to |p|c for a 
neutrino, the signal peaks at zero in C/ m i ss . 

The Pmiss distributions are shown in Fig. 3, where the tagged-P is required to be in the region of 
1.855 GeV/c 2 < Mbc <1-875 GeV/c 2 . DT yields, obtained by fitting the f/ m iss spectra, are listed in Table 3. 
Unbinned maximum likelihood fits are performed for all modes except for modes including cu. For modes 
including an w, binned least-square fits are performed to the w sideband-subtracted U m i ss distributions. 
In each fit, the Kev or Kyv signal is modeled by the MC-determined shape convoluted with a bifurcated 
Gaussian where all parameters are allowed to vary in the fit. Backgrounds for Kev are well described with 
a first-order polynomial. However, in the Kyv mode, backgrounds are more complex and consist of three 
parts. The primary background comes from D — > Ktttt 0 decay. To better control this background, we select 
a sample of D — > Ktttt 0 in data by requiring P ex tra >0.5 GeV, in which the t/ m iss shape of Kirn 0 is proved 
to be basically the same as that in the region of E ex tra <0.2 GeV in MC simulation. The selected Ktttt 0 
sample is used to extract the resolution differences in the t/ m iss shape of Ktttt 0 in MC and data, and to 
obtain the D Ktttt 0 yields in E ex tra >0.5 GeV region. Then, in fits to t/miss, the Ktttt 0 is described by 
the resolution-corrected shape from MC simulations and its size is fixed according to the relative simulated 
efficiencies of the E ex tra >0.5 GeV and P ex tra <0.2 GeV selection criteria. The second background from Kev 
events is modeled by a MC-determined shape. Its ratio to the signal yields is about 3.5% based on MC 
studies and is fixed in the fits. Background in the third category includes all other background processes, 
which are well described with a first-order polynomial. 

3. Systematic uncertainties 

Most sources of uncertainties for the ST or DT efficiencies, such as tracking, PID, and t r°, rj, Kg 
reconstruction, cancel out in determining yep- The main systematic uncertainties come from the background 
veto, modeling of the signals and backgrounds, fake tagged signals, and the CP- purity of ST events. 

The cosmic and Bhabha veto is applied only for the KK and tttt ST events which have only two tracks. 
The effect of this veto is estimated based on MC simulation. We compare the cases with and without this 
requirement and the resultant relative changes in ST efficiencies are about 0.3% for both the KK and tttt 
modes. The resulting systematic uncertainty on yep is 0.001. 

Peaking backgrounds are studied for different ST modes, especially for pir backgrounds in the Kg-rr 0 tag 
mode and iGy7r + 7r _ 7r° backgrounds in the K^uj tag mode. Based on a study of the inclusive MC samples, 
the fraction of peaking backgrounds in K^tt 0 is 0.3%. The uncertainties on yep caused by this is about 
0.001. Uncertainties from the sideband subtraction of peaking backgrounds for the Kgui mode are studied 
by changing the sideband and signal regions; changes in the efficiency-corrected yields are negligible. 

Fits to the Mbc and U m i ss spectra could induce systematic errors by the modeling of the signal and 
background shape. The MC-determined signal shapes convoluted with a Gaussian are found to describe 
the data well, and systematic uncertainties from the modeling of the signal are assumed to be negligible. 
To estimate uncertainties from modeling of backgrounds, different methods are considered. For the CP ST 
yields, we include an additional background component to account for the "0(3770) —> DD process with a 
shape determined by MC simulation whose yield is determined in the fit. The uncertainties in the fits to 
Mbc are uncorrelated among different tag modes, and the obtained change on yep is 0.001. For the DT 
semileptonic yields, the polynomial functions that are used to describe backgrounds in our nominal fits are 
replaced by a shape derived from MC simulations. For the Kyv mode, the size of the main background 
Ktttt 0 is fixed in our nominal fit, so the statistical uncertainties of the number of selected Ktttt 0 events 
introduces a systematic error. To estimate the associated uncertainty, we vary its size by ±1 standard 
deviation based on the selected Ktttt 0 samples. Systematic uncertainties due to the U m j ss fits are treated as 
positively correlated among different tag modes. We take the maximum change on the resultant yep , that 
is 0.006, as systematic uncertainty. 

The DT yields are obtained from the fit to the U m i ss spectra. However, one also has to consider events 
that peak at f7 m i ss but are backgrounds in the Mbc spectra, the so-called fake tagged signals. This issue is 


9 



Figure 3: Fit to the U m i ss distributions for selected DT events from data. In each plot, the solid line is the total fit, the dashed 
line in Rev shows the contribution of polynomial backgrounds, and the dash-dotted line in Kfiiy shows the contribution of the 
main Rmr 0 backgrounds. 
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Table 4: Summary of systematic uncertainties. Relative systematic uncertainties are listed for each tag mode in percent, while 
the resulting absolute uncertainties on y cp are shown in the last column. Negligible uncertainties are denoted by 



K+K~ 

+ 

K°ir°ir 0 

K°ir° 

K°u 


Vcp 

Background 

0.3 

0.3 

- 

0.3 



0.001 

Mpc fit 

0.4 

0.1 

2.4 

0.4 

0.1 

1.4 

0.001 

U miss fit (Kev) 

1.8 

1.3 

2.4 

1.6 

8.1 

1.2 

n nn£ 

Umiss fit 

3.2 

7.0 

4.6 

2.5 

1.7 

1.7 


Fake tag (Kev) 

0.2 

1.4 

0.9 

1.2 

3.1 

0.4 

n nn9 

Fake tag (Kyv) 

1.0 

0.7 

0.5 

0.9 

4.8 

0.4 


CP-purity 



0.4 


0.2 

0.2 

0.001 


examined by fitting to the Mbc versus U m i ss two-dimensional plots. From this study, the fake tagged signal 
component is proved to be very small. The resulting difference on yep is 0.002 and assigned as a systematic 
uncertainty. 

We study the CP-purities of ST modes by searching for same-CP DT signals in data. Assuming CP 
conservation in the charm sector, the same-CP process is prohibited, unless the studied CP modes are not 
pure or the initial C-odd D°D° system is diluted. The CP modes involving Kg are not pure due to the 
existence of small CPV in K° — K° mixing [6]. However, this small effect is negligible with our current 
sensitivity. Hence, Kgir 0 is assumed to be a clean CP mode, as its background level is very low. As a 
conservative treatment, we study DT yields of (Kgir 0 , Kgir 0 ) to verify its pure CP-odd eigenstate nature 
and the CP-odd environment of the D°D° pair. The observed numbers of this DT signals are quite small, 
and we estimate the dilution of the C-odd initial state to be less than 2% at 90% confidence level. This 
affects our measurement of yep by less than 0.001. The purity of the Kgir 0 mode is found to be larger than 
99%. Due to the complexity of the involved non-resonant and resonant processes in Kgir°ir° and Kgui, the 
CP-purities of these tag modes could be contaminated. We take the mode K + K~ as a clean CP-even tag 
to test Kgir°ir°, and take Kgir 0 to test KgOj and Kgr). The CP-purities of Kgir°ir 0 , Kgco and Kgy are 
estimated to be larger than 89.4%, 93.3% and 93.9%, respectively. Based on the obtained CP purities, the 
corresponding maximum effect on the determined yep is assigned as systematic uncertainty. 

Systematic uncertainties from different sources are assumed to be independent and are combined in 
quadrature to obtain the overall yep systematic uncertainties. The resultant total yep systematic uncer¬ 
tainties is 0.007. 


4. Results 


The branching ratios of K^e^v and K T yi^v are summed to obtain Bd CPt ->i = BDcp^^Kev+BDcp^^Kpv 
To combine results from different CP modes, the standard weighted least-square method is utilized [6]. The 
weighted semileptonic branching fraction Bd cp± ^i is determined by minimizing 


x 2 = Y ■ 


B 


Dcp±- 


>i-B 


Dcp ±^>1 , 


( a CP±) 


( 10 ) 


where a denotes different CP-tag modes and crg P± is the statistical error of £>g CP± _^ for the given tag 
mode. The branching fractions of Bd cp± ^,i and the obtained Bd cp± ^i are listed in Table 5. Finally, yep 
is calculated using Eq. (3), with Bd cp ±^i replaced by Bd cp ±^i- We obtain yep = (—2.0 ± 1.3 ± 0.7)%, 
where the first uncertainty is statistical and the second is systematic. 


5. Summary 

Using quantum-correlated D°D° pairs produced at s = 3.773 GeV, we employ a CP-tagging technique 
to obtain the yep parameter of D° — D° oscillations. Under the assumption of no direct CPV in the D 
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Table 5: Values of branching ratio of Dcp±—>i obtained from different tag modes and the combined branching ratio. The 
errors shown are statistical only. 


Tag Mode 

BDcp-^Kev (%) 

^Dc P— — yKfiP (%) 

B Dcp ^i (%) 

K+K~ 

3.44 ± 0.12 

3.33 ± 0.12 

6.77 ± 0.17 

i 

P 

+ 

1= 

3.29 ± 0.18 

3.35 ± 0.20 

6.64 ± 0.27 

KgTT°Tr 0 

3.40 ± 0.18 

3.48 ± 0.18 

6.89 ± 0.26 

^Dc p — —yl 

3.40 ± 0.09 

3.37 ± 0.09 

6.77 ± 0.12 

Tag Mode 

$Dcpp—tKeis (%) 

^Dcp-\-—^Kfiu (%) 

B Dcp+ ^ (%) 

K° s ir 0 

3.62 ± 0.10 

3.33 ± 0.10 

6.96 ± 0.15 


3.32 ± 0.21 

3.81 ± 0.21 

7.14 ± 0.30 

K%r, 

3.68 ± 0.26 

3.84 ± 0.29 

7.52 ± 0.40 

Bdcp+—>1 

3.58 ± 0.09 

3.46 ± 0.09 

7.04 ± 0.13 


sector, we obtain yep = (—2.0 ± 1.3(stat.) ± 0.7(syst.))%. This result is compatible with the previous mea¬ 
surements [18, 33, 34, 35] within about two standard deviations. However, the precision is still statistically 
limited and less precise than the current world average [6]. Future efforts using a global fit [36] may better 
exploit the BESIII data, leading to a more precise result. 
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